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ABSTRACT 



Supersymmetric model contributions to the neutron electric dipole moment arise via quark elec- 
tric dipole moment operators, whose matrix elements are usually calculated using the Naive 
Quark Model (NQM). However, experiments indicate that the NQM does not describe well the 
quark contributions Aq to the nucleon spin, and so may provide misleading estimates of electric 
dipole operator matrix elements. Taking the Aq from experiment, we indeed find consistently 
smaller estimates of the neutron electric dipole moment for given values of the supersymmet- 
ric model parameters. This weakens previous constraints on CP violation in supersymmetric 
models, which we exemplify analytically in the case where the lightest supersymmetric particle 
(LSP) is a U(l) gaugino B, and display numerically for other LSP candidates. 



The fact that the neutron electric dipole moment d n is expected to be so small in the 
Standard Model means that it is an excellent window for possible physics beyond the Standard 
Model, which may cast valuable light on the still-mysterious phenomenon of CP violation. Of 
special interest to us is the information that d n may provide on supersymmetry, which has been 
a focus of attention for quite some time [1], and in particular on the minimal supersymmetric 
extension of the Standard Model (MSSM). The MSSM contains many possible sources of CP 
violation beyond the Kobayashi-Maskawa phase of the Standard Model, whose presence could 
have significant impact on studies of supersymmetric candidates for cold dark matter [2] . 

Calculations of d n involve two aspects: a short-distance part that may be calculated per- 
turbatively in the Standard Model or one's chosen extension of it, and a long-distance non- 
perturbative QCD calculation of the matrix element (s) of the short-distance operator (s). In the 
case of the MSSM, it is generally thought that the value of (d n ) is dominated by quark electric 
dipole operators jQcr^qF^ [3], whose short-distance coefficient is related, in particular, to the 
chargino/neutralino sector that also features in dark-matter calculations. 

In the past, the matrix elements of these quark electric dipole moment operators have 
generally been estimated using the naive quark model (NQM). However, although this model 
is usually fairly reliable, there are indications that its predictions for the spin content of the 
nucleon are misleading [4]. In particular, there appears to be a non-zero strange quark contri- 
bution As to the nucleon spin, and the total quark contribution Au + Ad + As is quite small. 
These observations suggest that the NQM may not be adequate to describe the contributions 
of the quark electric dipole operators to the neutron electric dipole moment. In particular, 
a contribution from the strange quark could be significant, and there might be a systematic 
cancellation between the u, d and s contributions. This would mean that previous constraints 
on CP-violating parameters in the MSSM derived from upper limits on d n should be relaxed. 

Here we study this possibility analytically in the idealized case where the lightest super- 
symmetric particle (LSP) is a U(l) gaugino B, and numerically over the (//, M 2 ) plane that 
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characterizes the neutralino and chargino masses in the MSSM. We do indeed find a systematic 
suppression of d n compared to the NQM prediction for the same values of the CP-violating 
parameters in the MSSM. This in turn means that the corresponding upper bounds on these 
parameters inferred from the experimental upper limit on d n should be relaxed. This provides 
some breathing space for CP violation in the MSSM that may not be unwelcome, and in turn 
increases further the upper limit on the LSP mass inferred from calculations of its relic density 
that allow for CP-violating phases [2].* 

We first discuss the contribution of quarks to d n in terms of their contribution to the 
total neutron spin. The dipole operator ^qa^qF^ yields a contribution d q (Aq) n to d n , where 
( n l \qinl5q\ n ) = (Aq) n Sn defines the fractional contribution of the quark to the neutron spin 
S M , and d q is the fundamental quark electric dipole moment that arises due to CP violation. The 
fractions (Aq) n are related by isospin symmetry to the fractions Aq of the proton: (Au) n = Ad, 
(Ad) n = Au, and (As) n = As for the u, d and s quarks. In the non-relativistic NQM one 
would have Au = |, Ad = — |, and As = 0. Thus, d n = ^dj — ^d u in this case, which is the 
relation that has always been used in the previous literature on this subject. However, there are 
good indications that neither this model nor a simple relativistic modification of the NQM are 
adequate to describe the spin structure of the nucleon. Recent results [4] have confirmed the 
earlier result by the European Muon Collaboration (EMC) [5] that the total quark contribution 
Au + Ad + As is quite small. An analysis of the combined deep inelastic scattering and hyperon 
/?-decay data yields Au = 0.82 ± 0.03, Au = -0.44 ± 0.03, and As = -0.11 ± 0.03 [4]. We 
shall refer to these as EMC values of the fractions Aq. 

It is easy to see how the EMC values of the spin fractions tend to give smaller values of d n . 
Since d q oc m q , to a good approximation we have that d n oc m^Au + m s As. This is so because 
\Au\ ~ 2|Arf|, and ~ 2m u , therefore the u quark contribution is subdominant. Then, since 
the d and s quarks have the same quantum numbers, and m s m 150 MeV 3> m 10 MeV, we 

*Note that, in general, constraints based on the upper limit to the electron's electric dipole moment may be 
similar to those based on d n , if one assumes similar squark and slepton masses. 
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have contributions of approximately equal magnitude and opposite sign for small and negative 
values of As, as indicated by the EMC. 

We now turn to the calculation of d n in the MSSM to show that these expectations are 
borne out. The calculation of the fundamental quark electric-dipole moments in the MSSM 
has been extensively discussed in the literature. Here we shall follow the phase conventions 
and complete one-loop analysis of ref. [6] in order to calculate d n from the fundamental quark 
moments d q . 

We begin our discussion with the simple case in which the LSP is the pure gaugino B. This 
is a good approximation to the composition of the MSSM LSP in the region 3> M 2 of the 
parameter space characterized by the magnitude of the Higgs mixing mass \i and SU(2) gaugino 
mass M 2 (for simplicity, we shall assume GUT conditions on the gaugino masses so that we need 
to consider only one of them), and where the LSP makes an attractive dark-matter candidate. 
The quark contributions, d q , to d n = J2 q & n d g (Aq) n involve only gluino and chargino exchange 
if we neglect the possible contributions of heavier neutralinos. Note that in the MSSM there is 
no net contribution to d q from the heavier neutralinos in the region 3> M2, so this is not an 
unrealistic model. The contribution due to gluino exchange is given by (assuming a common 
squark mass; we shall relax this assumption in our later numerical analysis) 

where m q e tlq = R q /j, + A* is the off-diagonal term in the squark mass matrix, r = m~/M~, and 
K{r) = — [1 + 5r + 2r(2 + r) lnr/(l — r)]/2(l — r) 3 . Here nig(M q ) is the gluino (squark) mass, 
R q = cot/3 (tan/3) for an isospin T 3q = +^(— \) quark, and tan/3 is the usual ratio of the vacuum 
expectation values of the hypercharge +1 and hypercharge —1 Higgs scalars. We shall assume 
throughout a common, complex, dimensionless parameter A for the trilinear terms. Thus, if 
\A\ <C \n\ and sin^ ~ sin7^, 

d% (X J2Qg R q m q( A Q)n • (2) 
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In Figure 1 we plot, as a function of tan/3, the ratio of the values assuming NQM and 
EMC values of the spin fractions Aq. Note that the values are of opposite sign. As can 
be seen there, the ratio can be fairly large for tan/5 2 and for the central EMC value of 
As. Moreover, the ratio can be large for all values of tan/3 if As = —0.08, which is just la 
below the central EMC value: d n is fairly sensitive to the value of As because m s ^> mj. In 
the opposite case that \A\ 3> \n\ (or the off-diagonal terms m q e %lq are all similar if A is not 
common), oc Y,qenQq m q(^Q)n an d the ratio is large as well, \d^ N q M /d^ EMC \ m 4.3 (28 for 
As = -0.08). 

The chargino-exchange contribution to d q is given by (again, assuming a common squark 
mass) 

where 9^ is the phase of fi, Ti = rn^./M-, m W2 (m Wl ) is the heavy (light) chargino mass, q' 
refers to the doublet partner of q, I(r) = [1 + r + 2rlnr/(l — r)]/2(l — r) 2 , and J(r) = 
[3 — r + 2rlnr/(l — r)]/2(l — r) 2 . Since the factor in square brackets is of opposite sign and 
approximately equal magnitude for quarks of opposite isospin (differing by ^ 10% for r\ < 1), 
in this case we have approximately 

dC(xJ22T 3g R g m g (Aq) n . (4) 

In Figure 2 we plot, as a function of tan/3, the ratio of the values d% assuming NQM and EMC 
values of the spin fractions Aq. As can be seen there, the results are qualitatively the same as 
in Fig. 1. 

These results, for the simple neutralino sector we have assumed, help illustrate and explain 
the results of our more general numerical analysis of the MSSM, to which we now turn. There 
are four neutralinos that contribute to d q in this model in addition to the charginos and the 
gluino. We have diagonalized numerically the neutralino mass matrix to find the neutralino 
composition in the gaugino-Higgsino basis. From this composition, which now depends on 
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the parameters M 2 , tan/5 and 6 M , we determine the contribution of each neutralino to the 
quark moment d q . We do not assume a common squark mass as before. Rather, we use the 
squark compositions and masses determined by the diagonalization of the squark mass matrix. 
In addition, in order to have a more realistic squark mass spectrum, we allow for radiative 
corrections to the soft supersymmetry breaking masses assuming a common sfermion mass uiq 
at the GUT scale. Thus, we have differing left and right soft supersymmetry breaking masses. 
The ratio of d n calculated with NQM values of the spin fractions Aq to d n calculated with 
EMC values, \d^ M /d^ MC \, depends on many parameters now. In Figure 3 we exemplify our 
results in the usual (//, M 2 ) plane, for tan/3 = 2 and 9a = 90°. Here we would expect, on the 
basis of our previous results, to find a ratio ^ 1.7 throughout the plane. We show the ratio 
\dn QM / d n MC \ for A = m = 300 GeV with, in the top (bottom) panel, 0„ = 20° (160°). The 
light (darker, darkest) shaded areas correspond to regions where \d^ QM /d^ MC \ > 1.5 (3, 6). 
There are only light-shaded and unshaded areas in the region n < 500 GeV or M 2 > 500 GeV, 
and this is so throughout the (fi < 1000 GeV, M 2 < 1000 GeV) plane for much larger m . 

We find that only in a small fraction of the unshaded regions is the ratio smaller than unity. 
There is only a very narrow band in the unshaded region of the bottom panel of Fig. 3 (running 
parallel to the darkest shaded area) where the ratio is substantially smaller than unity. In that 
band, as well as in the darkest shaded areas of Fig. 3, there is a near cancellation between the 
chargino and gluino contributions. However, the near cancellation occurs over a much wider area 
for EMC values of the spin fractions because both contributions are made small, and therefore 
the difference of small numbers makes for a naturally small d n . For example, in the bottom 
panel of Fig. 3 \d% QM /d% MC \ < 0.1 only in a band of width Afi « 15 GeV for M 2 = 50 GeV, 
whereas \d% QM /d% MC \ > 10 in a band A/i « 80 (125, 750) GeV for As = -0.11 (-0.1, -0.08). 

Thus, our previous results are confirmed qualitatively by this more complete treatment 
of the MSSM. Moreover, we draw the reader's attention to the very dark areas, where even 
for tan/5 = 2 and As = —0.11 there can be a substantially smaller d n for relatively light 
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squark masses and large #4! This contrasts with the typical situation if one assumes the 
NQM of the nucleon structure. In that case, the low value of d n required by experiment 
< 1.1 x 10~ 25 e cm [7]) can only arise for relatively heavy (mass ~ TeV) squark masses if the 
CP-violating phases are 0(1) [6]. Such large masses, however, do not allow the neutralino LSP 
to be a viable dark matter candidate because its relic density would far exceed the observational 
constraints on the mass density of the universe. Recently [2], however, it has been pointed out 
that there is some room in the MSSM for a dark matter B even for relatively large CP-violating 
phases, provided B is sufficiently heavy. In Fig. 3, on the other hand, we see that even for a 
light neutralino (small M 2 ) there is the possibility of a sufficiently small d n for relatively light 
squark masses, and that d n can be small over a much wider region of parameter space. 

One possible source of uncertainty in the value of d n in the dark regions is that the other 
contributions might prevent d n from being small. However, on quite general grounds d n must 
be proportional to the neutron spin, and thus the only other quark contribution that might 
not depend on the spin fractions would depend on orbital angular momentum, which appears 
to be small [8]. 

We conclude on the basis of the above arguments that a more realistic model of the quark 
structure of the nucleon spin indicates that d n is smaller than previously thought in super- 
symmetric theories. This weakens constraints on CP-violating parameters in supersymmetric 
models, and further increases the upper limit on the LSP mass inferred from its relic density by 
calculations that allow for CP-violating phases. We have also shown that d n could be substan- 
tially smaller if the LSP is a relatively light B, thus raising the possibility that a light neutralino 
could constitute the cold dark matter without requiring particularly small CP-violating phases 
in the theory. 

The work of RF is supported by an NSF grant. He thanks the hospitality of the Theory 
Division at CERN while this work was started. 
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FIGURE CAPTIONS 



Figure 1 - The ratio —d^ N Q M /d^ EMC , obtained using Eq. (2), of the gluino contribution 
to d n in the simple model with a B LSP calculated with NQM values of the spin fractions 
Aq (Au = |, Ad = — |, As = 0) to the contribution calculated with EMC values (Au = 
0.82, Ad = -0.44, and As = -0.11 (-0.08) for the solid (dotted) curve). 

Figure 2 - The ratio —d^ N q M /d^ EMC , obtained using Eq. (4), of the chargino contribution 
to d n in the simple model with a B LSP calculated with NQM values of the spin fractions Aq to 
the contribution calculated with EMC values. The values of the spin fractions are as in Figure 
1. 

Figure 3 - The ratio \d% QM /d® MC \ of d n in the MSSM calculated with NQM values of the 
spin fractions Aq to d n calculated with EMC values. The values of the spin fractions are as 
in Figure 1. We have assumed 6^ = 20° (160°) in the top (bottom) panel. In the light (dark, 
darker) shaded areas, \d^ QM /d^ MC \ > 1.5 (3, 6). 
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FIGURE 3 



11 



